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In the claims: 

1 . A method for gamut mapping of an input image using a : 
algorithm, comprising: 

receiving the input 

converting color, 
electrical values set; 

applying the space 




fions of an image pixel set to produce a corresponding 
fg algorithm to the electrical values set to produce a 



color-mapped value set; and 

reconverting the color-mapped value set to an output image. 

Tim mi'lliiiil ol'i l.iini I win n in lln r |>nr r vniynip ,i I g o iitIi in mlnimiTn 
variational problem represented by: 

E(u)= \(d 2 +a|VD| 2 )/Q, subject to u e i9 , wherein Q is a support of the input 

n 

image, a is a non-negative real number, Dt^*(u-u 0 ) 5 g is a normalized Gaussian kernel 
with zero mean and a small varianc^^Uo is the input image, and u is the output image. 
3. The method of claim 2 5> tartJier comprising: 

solving the variational problem at a high value of a; 
solving the variational problem at a low value of a; and 
averaging th^ solutions. 

The methieKi of claim 3, wherein the step of averaging the solutions comprises 
using a spatially adaptive weighting scheme, comprising: 

U /nat [*> A = W [ k > j}* small [ k > 7 X 1 ~ W i k > J^high [ k > A 

/herein the weight w[k, j] 9 comprises: 
1 



4. 



w 



l + /?|V#* Wo 



-, and 



wherein 3 is a non-negative real number. 

The method of claim 2, wherein the variational problem is solved according to: 
du 



dt 



= ag* AD - g * Z>, subject to u e 3 . 



aim 2, wherein the space varying algorithm^ solved according 
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lr=g*g*( ( 4 n -u 0 ) 

7. The method ofp&im 1, wherein the space varying algorithm minimizes a 
variational problemrepresented by: 

E(u)=Mp l (D)+ap 2 (\VD\))dCl, subject to mg5, wherein p x and /? 2 are scalar 
n 

functions. 

8. TKe method of claim 2, further comprising: 
lecimating the input image to create one or more resolution layers, wherein the 

one o/ more resolution layers comprises an image pyramid; and 

solving the variational problem for each of the one or more resolution layers. 

9. / The method of claim 1, wherein the method is executed in a camera. 

— Ihe method o T claim 1, wherein the metl i ud is eAe cut e d4ft---gr^rinter. 

11. A method for color gamut mapping, comprising: 
converting first colorimetric values of an input image^ttTsecond colorimetric 

values of an output device, wherein output values jjpe"constrained within a gamut of the 
output device; and 

using a space varying^gomtinff that solves an image difference problem. 

12. A computer-readable' mjtemoiy For color gamut mapping, comprising an instruction 
set for executing color gamut mgppitfg steps, the steps, comprising: 

converting firsj^colorimetric values of an original image to second colorimetric 
values, wherein output values are constrained within a gamut of the output device;using a 
space varying algorithm that solves an image difference problem; and 

optimizing a solution to the image difference problem. 
13^__Jhe-^e nipuier-l ' eadable iiieuiuiy uf daimj g ^x^ image diff erence 

problem is 
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1 ^a bject to uk=3, Wh erein Li is a support 01 an input image, a is a n Oll mgalivi iml^ 

2 number, D=g*(u-Uo), g is a normalized Gaussian kernel with zero mean^ad^mall 

3 variance a, Uo is the input image, and u is an output image. 

4 14. The computer-readable memory of claim 12,^h6rein the instruction set further 

5 comprises steps for: 

6 solving the image difference problprfi at a high value of a; 

7 solving the image difference pj^blem at a low value of a; and 

8 averaging the solutions. 

9 15. The computer-readable niemory of claim 14, wherein averaging the solutions 
10 comprises using a spatially adaptive weighting scheme, comprising: 

« final fy, j] = 4*> fa small [ k > J X 1 " W [ k > 7 [*> A md 
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wherein the weight w[k, j], comprises: 
1 



w 



-,and 



l + J3\Vg*u 0 \ 

wherein p 'is a non-negative real number. 

16. The comp/uter-readable memory of claim 12, wherein the image difference 
problem is represented by: 

E(u) = lj(p x (D)+ ap 2 (\VD\))dQ i , wherein p } and p 2 are scalar functions, 
/n 

17. The computer-readable memory of claim 12, wherein the instruction set further 
comprises steps for: 

decimating the input image to create one or more resolution layers, wherein the 
one or more resolution layers comprise an image pyramid; and 

solving the image difference problem for each of the one or more resolution 

layers. 

18. TheJ computer-readable memory of claim 17, wherein the instruction set further^ 
comprises steps for: 

(aj initializingjijirst resolution^ 

Jalculating a gradient G for the resolution layer, the gradient Cj composing: 
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and a k =a 0 *2 2(k " 1) ; 



(c) calculating a normalized steepest d^ent value Lj= Lj.i-|a 0 *HNSD*G, wherein 
ja 0 is a constant; 

(d) projecting the value onto^nstraints Proja (Lj), wherein Proja (x)is a projection 
of x into a gamut 9; and 

(e) for a subsequent resolution layer, repeating steps (b) - (d). 

19. A method for image enhancement using gamut mapping, comprising: 
receiving a input image; 

from the inpyr image, constructing an image pyramid having a plurality of 
resolution layers- 
processing each resolution layer, wherein the processing includes completing a 
gradient iterajaon, by: 

calculating a gradient G; 
completing a gradient descent iteration; and 

projecting the completed gradient descent iteration onto contstraints; and 
omputing an output image using the processed resolution layers. 

20. / The method of claim 19, wherein the gradient G, comprises: 
G = A(u-u 0 )+a k (u-u 0 ) 9 

wherein u is the output image, Uo is the input image, and a is a non-negative real 
number. 

21/. The method of claim 19, wherein completing the gradient descent iteration 
comprises calculating: 

^ =SG ^(G*AG) + a A ZG 2 ) ;and 
L j = L j-i - Mo 'Mnsd' G > 
^herein /j nsd is a normalized steepest descent parameter, ju 0 is a constant, k is a number 

f resolution layers in the image pyramid, and j is d specific resolution layer. 
£2T^^ The method of claim 19, wherein proj earing the completed gradient descent 
iteration onto the constraints is given b> 
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*— Lj-Piojs tEjS ■ 

wherein Proj8(x) is a projection of x into a gamut S. 

23. The method of claim 19, wherein constructing thejmage pyramid, comprises: 
smoothing the input image with a Gaussianjeemel; 
decimating the input image; and 
setting initial conductive L 0 = max^Sp}, wherein Sp is an image with the coarsest 

resolution layer for the image pyra 

24. The method of claim 23, vtfherein the Gaussian kernel, comprises: 
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25. The method'of claim 19, wherein processing each resolution layer further 
comprises applying a space varying algorithm to minimize a variational problem 
represented byji 

i?(w)= j[p 2 + a VD\ 2 ]iQ, subject to u <e 3 , wherein Q is a support of the 

image, and/ f D=g * (u-Uo), wherein g is a normalized Gaussian kernel with zero mean and 
small variance a, Uo is the input image, u is the output image, and wherein a is a non- 
negative real number. 

26. /The method of claim 19, wherein processing each resolution layer comprises 
applyi/ig a space varying algorithm to minimize a variational problem represented by; 

E{u)= j(p l (D) + ap 2 (\VD\))dCl 9 subject to «e5, wherein p x and p 2 

are Scalar functions. 

27. / The method of claim 26, wherein p l and p 2 M^ chosen from the group 



22 comprising p{x) -■^|_and^p^)= Vl -Kx^ 
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